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ABSTRACT 

An algorithm to model the effects of internal motion by the 
solar array and the high-gain antenna on the attitude of the 
Landsat-D spacecraft is presented here. The relative torque 
and angular momenta arising from the internal motions are 
assumed to be attitude- independent but are considered to be 
a source of attitude perturbations. The equation of motion 
for the three-body problem is derived and then compared with 
the one-body case. The effect of the internal motion on the 
control of the spacecraft is shown in a computer study of 
the problem. 

1 . INTRODUCTION 

The paper presents algorithms for modeling the effects of 
internally moving parts on the attitude of the Landsat-D 
(LSD) spacecraft. The internal motions considered here in- 
clude the rotations of the solar array to follow the Sun and 
the gimballed high-gain antenna to communicate with the 
Tracking and Data Relay Satellite (TDRS) (Reference 1) . The 
LSD system is treated as a rigid three-body system for de- 
scribing the equation of motion. Modeling the disturbance 
torques produced by moving appendages is very important for 
missions such as Landsat-D, which require accurate knowledge 
of the attitude and precise control of the spacecraft. 

The relative torques and angular momenta arising from the 
internal motions are considered as attitude-independent 
variables and as a source of attitude perturbations. The 
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external disturbance torques and the angular momenta caused 
by the internal motions are generated in a profile program 
(called PROFILE) on an IBM S/360-95 computer, where null 
attitudes are assumed and are transmitted to a truth model 
on a DEC PDP-11/70 computer that simulates the effects on 
the attitude. 

In this discussion, nonstandard rotations such as a 45-de- 
gree slew of the solar array to avoid interference with the 
antenna and the switching motion of the antenna from one 
TDRS to another are neglected. In addition to the rota- 
tional motions of the solar array and the antenna, the LSD 
spacecraft contains moving parts such as the thematic mapper 
and multispectral scanner (Reference 2). However, these 
motions are disregarded here because the motions are oscil- 
latory with a high frequency (=7 Hertz) and because they 
generate zero average angular momenta. 

Section 2 discusses the mathematical derivations of the 
equation of motion and pertinent terms such as the moment of 
inertia (MOI) tensor and the center of mass (CM) . When pos- 
sible, these terms are compared with the form for the one- 
body system used by the Multimission Modular Spacecraft 
(MMS)/Solar Maximum Mission (SMM) spacecraft. Section 3 
provides simulation results to compare the three-body and 
one-body cases. Conclusions resulting from the study are 
presented in Section 4. 

2. ANALYTICAL CONSIDERATIONS 

This section presents the mathematical modeling to describe 
the dynamic effects of the moving parts on the motion of the 
spacecraft. The equation of motion for the LSD mission is 
referenced at the CM of the entire system but is represented 
in a coordinate system that is fixed in the main vehicle. 

The CM of the entire system is calculated as a function of 
time. The MOI tensors for the moving parts are reevaluated 
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with respect to a set of time-independent axes parallel to a 
set in the main vehicle. Also calculated are the angular 
velocity of the appendages and the perturbation in the ex- 
ternal torques due to the changing positions of the append- 
ages. A comparison with the one-body problem is made. 

2.1 COORDINATE SYSTEMS AND TRANSFORMATIONS MATRICES 

The system under consideration, shown in Figure 1, consists 
of the main carrier vehicle, designated as body B^, and 
n(=2) moving bodies Bj (j=l, n) . Several coordinate sys- 
tems are convenient for discussing the relative motions. 
These are as follows: 

• Geocentric Inertial Coordinate System (GCI) (Refer- 
ence 3) 

• Orbit-Defined Coordinate System (OCS) where X 
(roll) is nearly along the spacecraft velocity vec- 
tor , Y (pitch) is along the orbit normal vector, 
and Z (yaw) is along the nadir vector 

• Spacecraft-Fixed Coordinate System (BCS) , which is 
fixed in the main vehicle B 

o 

• Coordinate systems fixed in moving parts such as in 
the solar array (SACS) or in the high-gain antenna 
(ANTCS) 
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NOTE: 



the CM of Bj 

the CM of the entire system 

the CM of Bj from r^j^ 
the hinge point of B^ 

the angular velocity of B^ in inertial space 

the angular velocity of B • relative to the main 

body Bq (oTj = 11^0 + 'uTj) 


Figure 1. Partitioning of the Satellite Into Main Body 
and Moving Parts 
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The transformation matrices (TRMA) to be used in this paper 
are defined as follows; 

1. TRMA from GCI to PCS ; [0] 


X 


l) X 



X 




X 

1 

<> 



X 




J 


where and denote the spacecraft position relative to 
the Earth and velocity unit vectors in the GCI frame, re- 
spectively. 


2 . Attitude direction cosine matrix from the PCS to 
the BCS : [A] . In the PROFILE Program [A] is given by the 

identity matrix because null attitudes are assumed. In the 
truth model, it is represented as 



( 2 - 2 ) 


using the small angle approximation, which is sufficient and 
valid, since only small perturbations are assumed; r, p, and 
y denote roll, pitch, and yaw angles in radian units, re- 
spectively. 
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(SA) 

3. TRMA from BCS to SACS : [C ’] . The solar array 

rotates around the y-axis and is driven to follow the Sun. 
Thus, its orientation is determined from the Sunline angle a 


[C 


(SA) j 



cos a 
0 

sin a 


0 -sin a 

1 0 

0 cos a 


(2-3) 


/N T . , 

Given the Sun unit vector, S = (S , S , S ) , in the BCS, 

X y z 

the rotation angle a is given by 



(2-4) 


because the Sun vector is perpendicular to the x-axis of the 
SACS . 

(ANT) 

4. TRMA from BCS to ANTCS ; . [C ] . The antenna 

has two gimbals with the inner gimbal angle, g^, repre- 
senting the elevation angle and the outer gimbal angle, 
g^, representing the azimuth angle. The orientation of 
the antenna is determined from the gimbal angles 


^11 

y L J z 


cos g^ cos sin g^ cos q^ -sin g 2 


(2-5) 


= -sin g^ cos g^ 0 

cos g^ sin g 2 sin g^ sin g 2 cos g 2 

The unit vector pointing from the spacecraft to TDRS is re- 
presented by P where P = (P„, P,,, P„)"^ in the BCS. 

X y Z 


14-6 



The gimbal angles are thus given by 


and 



= tan 


(P /P ) 

y/ jj/ 



(2-6a) 


(2-6b) 


Since g^, should align the antenna boresight (the x-axis 
in ANTCS) with the normalized pointing vector (P can be 

obtained from the spacecraft and TDRS ephemer ides . ) 

2 . 2 ANGULAR VELOCITY OF MOVING PARTS 

The angular velocity of the moving parts is used to calcu- 
late the internal angular momentum of the spacecraft for use 
in the equation of motion. It is easily seen from Equa- 
tion (2-3) that the angular velocity of the solar array is 
as follows: 




(2-7a) 


The time derivative of the rotation angle a can be com- 
puted numerically 


dg _ g(t) - g(t - At) 
dt At 


(2-7b) 
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Using Equation (2-5) the angular velocity of the high-gain 
antenna is 



(2-8a) 


where 


(t) - g^ (t - At) 

dt ” ~St 


(2-8b) 


For SMM, the angular velocity of the moving parts was not 
calculated . 


2.3 CENTER OF MASS 

For LSD, the CM of appendage in the BCS is given by 


Qj (t) 



Xj) 


H. X, 


(2-9) 


where Q^q represents the CM of B^ at the initial time (see 
Figure 1) . The rotation (or hinge) point is denoted by 

^ ^ J 

and QjQ - Xj represents the CM of B^ from the hinge point 
at the initial time. Then, at any later time, the CM will 
be represented by the first term of the right-hand side of 
Equation (2-9). The CM of each appendage changes as a 
function of time because the high-gain antenna rotates to 
track the TDRS, and the solar array rotates to track the 
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Sun. Consequently, the CM of the system, changes in 


time and is represented 

by 




n 

/ 





(2-lOa) 


r=0 / 

O 

II 



and the position of the CM of each appendage with respect to 
the CM of the system is 

(t) = "Q. (t) - (2~10b) 


For SMM, the CM of the system was fixed in time in the BCS. 

2.4 MOMENT OF INERTIA TENSOR OF THE SYSTEM 

The MOI of the system, [I^] , relative to axes parallel to 
the BCS axes passing through r^^^ is expressed by 


= S /a-r }l'3r<'=> * «ln, 

r=0 ' 

- [^^(t) + p^]^j 


( 2 - 11 ) 


where is the position vector of the mass dm^ of body 
relative to the CM of and the subscripts 1 and m 
represent the 1 and m components of the vector or tensor. 
Note that because is time-dependent, also de- 

pendent on time; in the remainder of this paper, the ex- 
plicit time-dependence will be dropped. 
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The above equation can be written as 


= Z) K - <’r' l '9r> 


r=0 


ra 




Imj 


( 2 - 12 ) 


since 


dm p = 0 
r 


~ ( r ) 

[I ] is the MOI tensor of represented in the BCS frame 
but relative to the CM of B^: 




(2-13) 


( T \ 

where [!' '] is the MOI of B^ represented in the coordinate 
system fixed in B^, . Equation ’(2-12) can be simply reexpressed 
by 


•V = i: 

r=0 




(2-14) 


with 




(r)' 


Im 


( 2 
+ M q ' 


<3 5 , - (q ) 1 (q ) ( (2-15) 

r p^r Im ^r 1 ^r m ) 


For the one-body problem, as represented by SMM, I is de- 
fined to be a constant in time. 

2.5 EXTERNAL TORQUES 

TWO external torques are discussed: the gravity gradient 

torque and the aerodynamic torque. The solar radiation 
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torque is similar to the aerodynamic torque, and the other 
external torques are not sensitive to the three-body problem. 

The gravity gradient torque, can be computed by 


^GG ~ 


n 


r = 0 


R + p 4- q 

(q; + rj X dm^ (2-16) 


where y is the Earth gravitational constant (“3. 986005 x 
14 3 2 -^ 

10 m /sec ). R is the spacecraft position vector from the 

Earth. Considering that 1^1 >> 1^ + ^ I, is simply, 

L r CjG 


N 


GG 




P,) 


r=0 

n 


3u V '' 1 ^ -Ci ^ r~(r), ^ 

= -f 2 _^ q^ X R(q^ • R) + R x [I ' M R 




n 


In 

"} 

R' 


R X [J R 


r = 0 


= ^ (R X [I ] R 
R-^ 


2-17) 


The expression for the one-body system has the ‘same form 
except for the replacement of [I^] by the constant [I]. 


To simplify the calculation of the solar radiation and aero- 
dynamic torques, the LSD spacecraft is modeled as an as- 
sembly of a cylinder for the main vehicle, flat plates for 
the solar array panels,, and a sphere for the antenna. Only 
the aerodynamic torque is discussed here because the modi- 
fications to the center of pressure (CP) are common in solar 
radiation and aerodynamic torques. 
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The aerodynamic torque, is 



(2-18) 


Here, "v denotes the spacecraft velocity unit vector, n^ 
denotes the normal unit vector for the ith surface, "q . 

f ± 

denotes the CP of the ith surface from ‘r’^^, p denotes 
the atmospheric density, and Cj^ denotes the drag coeffi- 
cient. The normal vectors, n^, for the solar array and 
antenna surfaces are dependent on time by 


/S 

n . 


(i) 




ft. 


lO 


(2-19) 


where n. 

lO 

sur face . 
puted by 


represents the initial normal vector for the ith 

'q . for the solar array and antenna are com- 
^cp , 1 


with 


cp,i 


= Q 


cp,i 


’CM 


( 2 - 20 ) 


cp,i 


-[< 


(i) 


(Q. 


cp , lO 


- ^i^ 


+ X. 


( 2 - 21 ) 


More consideration is required to specify for the 

solar array surfaces that are canted. The transformation 
matrix from BCS to these surfaces, is given by 


( 2 - 22 ) 

with the canted angle 9 . 
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For 

2.6 

The 

the 


the one-body case of SMM, n. and ^ are constants 

X cp,i 


EQUATION OF MOTION 

equation of motion for the LSD spacecraft is written 
form 


in 


IX 

dt 


= f (Y(t) , 


t) 


(2-23) 


where Y = (q^, (M = If 2 , 3, 4) denotes the 

Euler symmetric parameters representing a rotation from the 

GCI to the spacecraft-fixed coordinate frame, is the 

total angular momentum of the spacecraft, and is the 

w 

wheel momentum. 

The body angular momentum of the main vehicle, is 

B 

given by the total spacecraft angular momentum minus the sum 
of the wheel momentum, payload momentum, IlL , and the angu- 
lar momentum, caused by the internal motions 



(2-24) 


Lg depends on the angular velocity of the main vehicle, 
Ioq, and depends on the angular velocity of moving 

parts, . To formulate these mathematically, the 
angular momentum of the total system, L^, ignoring wheel 
and payload momenta, is considered 


^T S ^INT 


n 


= ^ / (<3^+ X (^j. X ^^) dm^ 

r=0 J 

n 

= ^ } Mj. X q'j. + (% + 1^;.) 

r=0 ' 


(2-25) 
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with some computation, can be shown as 




= [I 


S + CK^'^h oj: 


r=0 


where 


[K^ hin^ - Mr jqj. • - X^) - (r^j^ - X^,)^ 


Thus, the body rate of the main carrier is simply 


o)o = 


and caused by the internal motion, is 

INT 


'INT 


r = l 


The time derivatives of the Euler symmetric parameter, 
q , can be obtained as 

dq , 

dt^ = 2 ^yv 


with 




1 

o 

3 

"“2 

r 

3 

-a >3 

0 

“l 

“2 

"^2 

1 — i 

3 

1 

0 

<^3 

1 — \ 
3 

1 

J 

CM 

3 

1 

-(03 

1 

3 


(2-26) 


(2-27) 


( 2-28) 


(2-29) 


(2-30) 
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The time derivative of the total angular momentum of the 
spacecraft is given by the Euler equation as 

For SMM, the body angular momenta, 'L'g, is given by 




( 2 - 32 ) 


with the payload momentum L„. The spacecraft body rate, 

K 

1 ^, is determined by 



where [I]~^ is the inverse of the spacecraft MOI tensor. 
The time derivatives of the Euler symmetric parameters, 
q , can be obtained as 


dq 

dt 


1 

2 


' )JV 



( 2 - 34 ) 


The time derivatives of the total angular momentum of the 
spacecraft are given by the Euler equation as 

3^ ' Next * It X “ ■ <2-33) 

with the external torque, 
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3. SIMULATION RESULTS 


A computer study of the effect of the three-body problem on 
the motion of the spacecraft has been made using the general 
equations derived here. Since the spacecraft is subject to 
noticeable external torques, a control law that provides 
compensatory torques was necessary to keep the spacecraft 
near null attitude. The one-body case, using the same con- 
trol law, was also studied. 

The roll, pitch, and yaw of the spacecraft main carrier for 
both cases is shown in Figures 2 through 4. The results of 
the three-body case are represented by the "X” points and 
the results of the one-body case are shown as open circles. 
Note that both cases are subject to the same control law. 
This control law attempts to make the pitch, roll, and yaw 
zero and to bring the spacecraft rate to null. This control 
law is the same one (Reference 4) that Landsat-D will use 
during its acquisition phases. The torque applied to each 
reaction wheel is as follows: 

for the roll axis. 


T = K (k Ar + 0 ) ) (3-la) 

r r r r 

for the pitch axis, 

Tp = Kp [kp(Ap + B) + tOp] (3-lb) 

and for the yaw axis, 

where Ar and Ap are the roll and pitch attitude errors as 
determined by an Earth sensor; K , K , K. , k , k , k , and k 

L ^ y L y 
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Figure 2. Spacecraft Roll Versus Time 
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Figure 3. Spacecraft Pitch Versus Time 
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Figure 4. Spacecraft Yaw Versus Time 
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are constants; B is a bias to compensate for the orbital ro- 
tation; and (ji , ui , and co are the angular velocity along 
the roll, pitch, and yaw data. Because of the values used for 

k , k , and k , the control law is much more sensitive to 
r' p' y' 

the spacecraft rate than to the attitude error. 

Most of the structure seen in the plots is a result of the 
control law. However, since the control law is the same, 
the differences in the plots are a result of the three-body 
problem. Note in Figure 2 that after 4.5 minutes the con- 
trol law has the roll rate to zero for the one-body problem 
but not the three-body problem. Likewise, after 2.5 min- 
utes, the pitch rate of the one-body problem is under con- 
trol. 

4. CONCLUSIONS 

The conversion of the rigid one-body problem to the three- 
body problem has added another dimension to the study of 
dynamics. Although the exact perturbations in motion are 
obscured by the control law used, the effects are still im- 
portant in control of the spacecraft. 

The algorithms used in this paper can be applied to other 
spacecraft such as the Space Telescope to study important 
low-frequency effects, as in this paper, and also higher 
frequency effects that will cause jitter. 
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